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Abstract. 

KTa03-based solid solutions exhibit a variety of interesting physical phenomena. 
To better understand these phenomena, we performed first-principles calculations on 
[Ki.^LijTaOa (KLT) and K[Tai_ x ,Nb x ]C>3 (KTN) supercells. Our results show Li 
displacements and potential barrier heights in KLT that are in excellent agreement with 
values obtained from experimental fits. Dramatic changes in B-site dynamical charges 
occur in KTN in response to changes in near neighbor (nn) coordination. These effects 
can be explained by heterogeneity in the local electronic dielectric permittivity. 



INTRODUCTION 

KTaC>3 (KT) is one of the rare perovskite-type materials that remains cubic 
down to zero Kelvin. It is also interesting because it is a quantum paraelectric 
and, therefore, its ferroelectric phase transition at low temperatures is suppressed 
by zero-point quantum vibrations. KNbC>3 (KN) is ferroelectric even at room tem- 
perature and exhibits a series of ferroelectric phase transitions with decreasing 
temperature. A variety of interesting phenomena occur when KT or/and KN form 
solid solutions via isovalent substitution of A or B cations; [Ki_ x ,LijTa03 (KLT) 
and K[Tai_ x ,Nb x ]03 (KTN) are solid solutions which exhibit both order- disorder 
and soft-mode characteristics. In this work, we use first-principles supercell cal- 
culations to investigate lattice dynamics, ground-state structures, and dynamical 
effective charges in KLT and KTN. 

Recent developments in quantum theory of polarization (see the review in [1]) al- 
low one to make first-principle (FP) computations of polarization and its derivative, 
dynamical charge. Such computations show that dynamical charges are anoma- 
lously large in ferroelectric perovskite-type crystals (see e.g. Ref. [2]), and this has 
a significant effect on the soft-mode behaviors of these materials. Variations of 
dynamical charges as functions of local chemical coordination, as occurs in solid 
solutions, have not been extensively studied. We present systematic results on 



the dynamical charges as functions of chemical configurations in KTN. The same 
tendencies were found for PMN supercells. 

METHODS 

We use the Vienna ah initio simulation package (VASP) [3,4], which treats elec- 
tronic structure within the framework of density function theory. Detailed compu- 
tations were performed for 40 ion supercells and a limited number of computations 
were done for an 80 ion supercell. The basis vectors were [002], [020] and [200] for 
the 40 atom cell (KLT40), and [220], [022] and [202] for the 80 atom cell. The k- 
vector mesh was constructed from 4x4x4 Monkhorst-Pack grid. Verification of this 
computational scheme was performed for pure KTaOs; the self-consistent lattice 
constant a = 3.96 A agrees with previous LDA and GGA computations [5]. The 
calculated lattice constant is smaller than the experimental value a = 3.983 A [6]. 

VASP computes interatomic forces and total energies for crystals, and allows 
global relaxation as well as constrained relaxation of chosen internal coordinates 
and/or lattice parameters. We used frozen-phonon methods to obtain force con- 
stants for computing lattice dynamics. To compute dynamical charges, we used 
Berry phase analyses as implemented in the VASP code by Martijn Marsman. 

KLT RESULTS 

We calculated the force constants matrix in KLT40 containing one Li per cell. 
The reference structure for the force constant calculations was obtained by first 
placing all ions on ideal perovskite positions and then optimizing lattice parameters 
and internal coordinates while maintaining full cubic symmetry. The resultant cell 
had a = 2 x 3.956 A. 

From the force constants, we obtained the dynamical matrix, and the normal 
mode frequencies and eigenvectors. Symmetry analysis of KLT40 shows 14 
triplets, one which is acoustic and 13 which represent polar transverse optical (TO) 
modes. In Table 1, we give the TO normal mode frequencies. 

TABLE 1. Computed TO normal mode frequencies v (in cm -1 ) for cubic 
K 7 LiTa 8 C>24 (KLT). Also shown are relative amplitudes of each mode in 
the fully relaxed tetragonal KLT structure a re i and relative contributions 
of each mode to the tetragonal phase polarization P re i 
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For comparison, we also calculated the force constants of a 40-atom cell of pure 
KT at the same lattice parameter. The TO normal mode results are shown in 
Table 2. The force constants in the two cases are nearly identical. The only inter- 
atomic force constants which change by more than 0.41 eV/A 2 are those involving 
the 12 O ions nearest the Li site, and which favor Li moving against its nearest 
neighor O sites. In fact, the one lattice instability found in KLT is essentially Li 
motion (85% of the dynamical matrix eigenvector) opposite its nearest neighbor 
O (12% of the dynamical matrix eigenvector). Because KLT40 and KT have such 
similar force constant matrices, it is not surprising that the normal mode spectra 
are very similar: note that each LO mode in KT has a corresponding KLT40 mode 
with very similar frequency. Thus, the FP results confirm a picture of KLT in which 
Li are displaced off-center. We emphasize, however, that a quantitative description 
of the energetics and dynamics of KLT requires that one also take into account the 
coordinated motion of the ions surrounding the Li. 



TABLE 2. Computed TO normal mode fre- 
quencies v (in cm -1 ) and dynamical matrix 
eigenvectors for KTaC-3 at a = 3.956 A. 
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The largest reduction in energy in KLT40 occurs when Li is displaced along a 
[001] vector. Our results thus support a picture of KLT in which each Li is displaced 
in one of the 6 [001]-type directions. We performed a full ionic relaxation of KLT, 
polarized along z. A breakdown of the displacement pattern into normal mode 
coordinates is shown in Table 1. Although there is some mixture of higher frequency 
modes due to anharmonic coupling, the displacement pattern is dominated by the 
eigenvector of the lattice instability. A more interesting picture emerges when 
we investigate the contribution of each mode to the polarization of the relaxed 
structure. Because the 121 cm" 1 mode corresponding to the soft mode of pure KT 
has such a large mode effective charge, it contributes 23% of the total polarization 
of KLT40. 

We investigated Li energy as a function of displacement direction and found that 
the smallest barrier between [001] wells are saddle points in the [110] directions. 
Fig. 1 shows the energy as a function of Li displacement d in the [001] and [110] 
directions (see also [7]). Since the barrier energy is comparable with the thermal 
energy, Li can jump between neighboring [0,0,3] wells via [110] saddle points. 

If all ions except for the Li are fixed, the Li displacements and potential barriers 
are very small relative those deduced experimentally (by fitting a mean field expan- 
sion to experimental data [8]) (upper curves in Fig. 1, see also [7]). However, fully 
relaxed calculations (lower curves in Fig. 1) yield good agreement: the equilibrium 
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FIGURE 1. Li adiabatic and single-ion potential relief in KTaC>3 in [001], [110] directions (a) 
and in the direction between the Li wells (b) 



Li displacement is [0,0,0.99 A] with a = 3.96 A and [0,0,1.07 A] with a = 4.0 A; 
the corresponding adiabatic potential barriers are ~103 meV (1190 K) and ~1 17 
meV (1360 K), respectively. The calculated values are somewhat larger than the 
experimentally estimated barrier of 86 meV (998 K) [8] , but they agree rather well 
with the 1200 K barrier reported by [9]. Computed barrier heights may be system- 
atically larger than fitted values if the latter are reduced by Li-tunneling between 
neighboring wells over excited states. For the 80 atom supercell we obtained the Li 
displacement 1.01 A and the energy barrier U =118 meV (1369 K), which is also 
somewhat larger than the experimental value of 86 meV (998 K) [8]. 

Fig. 1(b) represents our computations of Li potential relief along a vector con- 
necting neighboring Li-wells, e.g. [0,0, d] and [0,d,0] via a minimum energy [011] 
saddle. We computed the "frozen coordinate" (i.e all ions frozen except for Li) en- 
ergies, as well as the adiabatic energies (i.e with relaxed coordinates of the non-Li 
ions). Clearly, ionic relaxation promotes Li- hopping by reducing potential barriers. 

In a relaxed supercell, cations displace in the same direction as Li, and anions 
displace in the opposite direction. This result supports the idea of coupling be- 
tween Li displacements and ferroelectrically active optical host ion modes [10]. 
The largest local deformation is exhibited by the four oxygens that are closest to 
Li. For example, in the 80 atom cell where the Li-displacement is [0,0,1.01 A], 
their displacements are [0,-0.074 A,-0.074 A]. In the 40 ion supercell these dis- 



placements are similar, [0,-0.093 A,-0.093 A]. Such distortion was proposed [11] 
to explain experimental data on the temperature dependence of photocurrent in 
[Ki_ x ,Li x ]Ta03 solid solutions. 

These results are in excellent agreement with computations for LiTaC>3 and 
LiNb0 3 [12], which also indicated a large oxygen displacements similar to those in 
our calculations. These oxygen displacments were explained in terms of a reduction 
in elastic energy that occurs because oxygens prefer to rotate about B site cations 
rather than directly approaching them [12]. 

Based on the fact that Li hop over a saddle point, we derived the following 
expression describing the temperature dependence of the relaxation time: 

r = A J e( u+ay2 V kBT dy = aJ —^—e u( - T ^ /kBT = r (T) e W (1) 

— oo 

where: r is the relaxation time; ks is the Boltzmann constant; a is a constant. 
The temperature dependence of the prefactor is a consequence of the existense 
of different paths when Li overcomes the potential barrier. The temperature de- 
pendence of the average barrier, U(T), is governed by the dynamics of the ions 
surrounding a Li. This dependence can be understood if one introduces spatial and 
thermal potential barrier fluctuations 

e -U(T)/k B T = g£ 7 dx7 dv e (Uo+x+y)/k B T e -a 2 x 2 -P 2 (y/k B T) 2 

nk B T J J il /r)\ 

— oo — oo y£ ) 
— e -U /k B T+l/Al3 2 +l/4a 2 k 2 B T 2 

From comparison of the left and right sides of this equality one finds U (T) = 
Uq — k B T/Al3 2 — l/4a 2 /c B T. At high temperature the potential barrier decreases 
linearly with temperature because of thermal fluctuations (in excellent agreement 
with experiments performed for [K .957,Li .043]TaO3 [13]). At low temperatures 
the barrier height decreases because of spatial fluctuations. Also, critical slowing 
down of relaxation is common, at low temperature, but in KLT, Li-tunneling is a 
competitive process [14]. 

We also performed computations for 80 atom supercells with Li-Li pairs (Fig. 2). 
A Li-Li pair was introduced into an 80-ion supercell such that nn A-sites (separa- 
tion vector az) were occupied by Li. We considered the following Li-displacement 
configurations for this pair: (zz), (xx) and (xx). The energy required to convert 
a (zz) to (xx) configuration is 0.27 eV (2990 K) in excellent agreement with 
the reported barrier height of 2800 K [9,13], attributed to the rearrangement of Li 
pairs [15]. Significantly, the excitation energy required to move Li pairs into the 
(xx) configuration, 0.26 eV, is close to the energy required to excite the (xx) con- 
figuration. This occurs because of the indirect dipole-dipole interaction between Li 
ions in Li-Li pairs over the KT-like soft mode. 

The computed dynamical charge of Li in KTaC^Li is 1.2 e at the equilibrium 
position but it becomes lower than 1.0 at small Li displacements. 
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FIGURE 2. A-site configurations in a [K 1 4/ 16 Li 2 / 16 ]Ta03 supercell (undecorated sites are oc- 
cupied by K): a) (zz), b) (xx), c) (xx) 



KT AND KN RESULTS 



The dynamical charges of Ta (Z^ a ) and Nb (Z^), in pure KT and KN, respec- 
tively, were calculated as functions of [001] displacement in cubic structures with a 
= 4.015 A, using a 6x6x6 Monkhorst-Pack k-point grid. The results are shown in 
Fig. 3. 

We also computed electrostatic potentials on the ions and their derivatives. Dis- 
placement induced reductions in the second derivative of the potential were con- 
sistent with changes in Z^ a and Z* Nb . Indeed one can write a general relation: 
Ei{u) = £,i(u)Z*Ui where ^j(O) ~ l/^oo and £ ro is a high frequency dielectric con- 
stant. Hence the reduction in local field (and charge-transfer [1]) that is caused by 
ionic displacement is responsible for the calculated decrease in dynamical charge. 
Note that this decrease is rather large: e.g. for Nb in KN Z* m = 9.8 at u=0; 8.5 at 
u=0.1 A; and 6.8 at u=0.2 A. 



KTN RESULTS 

Results of Z* computations for KTN supercells are listed in Table 3 and some 
results for [001] ordered supercells are shown in Fig. 4. Calculations for different 
[hkl] ordered superstructures, (where [hkl] indicates an ordered sequence of layers 
perpendicular to the [hkl] vector) were all done with an 8x8x8 Monkhorst-Pack 
grid. All computations were performed for the cubic structure with the fixed lattice 
parameter a = 3.96 A. This was done in order to show the pure effect of the 
neighbors on the dynamical charge. The effect of the ionic displacements were 
discussed in the previous section. 

In the cubic KN, Z* Nb ^ zz = 9.9, which is larger than Z^ azz = 8.8 in KT. Our 
results are consistent with the studies of Singh [16] which compared covalencies in 
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FIGURE 3. Dependence of the dynamical charge on the cation displacement: 1. KN, 
a = 4.015A. 2. KT a = 4.015A. 3. KT a = 3.98A. 



TABLE 3. Dynamical charges in KTN supercells: the first two columns represent the 
relative numbers of the Nb, x^jb, and Ta, XTa, planes; the vector [hkl] shows the direction 



perpendicular to the planes. 



X7V6 


XTa 


[hkl] 


Z Nb.zz 


Z* Nb,xx 


Z*Ta,zz 


& Ta,xx 


1 







9.91 


9.91 









1 








8.79 


8.79 


1/2 


1/2 


[001] 


9.06 


9.88 


9.56 


8.81 


1/3 


2/3 


[001] 


8.78 


9.88 


9.30 


8.80 












9.31 


8.81 


1/4 


3/4 


[001] 


8.62 


9.88 


9.18 


8.80 












9.23 


8.80 












9.16 


8.80 


2/3 


1/3 


[001] 


9.31 


9.89 


9.89 


8.82 








9.31 


9.89 






3/4 


1/4 


[001] 


9.47 


9.89 


10.06 


8.81 








9.41 


9.90 












9.47 


9.89 






1/2 


1/2 


[111] 


8.75 


8.75 


9.22 


9.22 


1/2 


1/2 


[110] 


9.87 


9.23 


8.84 


9.68 




1 2 3 1 2 3 



Number of Nb planes Number of Ta planes 

FIGURE 4. Dependence of the dynamical charge on the number of Ta- and Nb- planes in KTN 
(001) supercells 



KT and KN. In layered Nb/Ta superstructures Z* m ^ zz decreases with the number 
of nn Ta-planes, and Z^ azz increases with the number of nn Nb-planes. Both 
Z* B zz s approach saturation values; Z* B ^'s, however, exhibit no significant changes. 
These results are similar to those of previous computations for Nb/Al supercells 
in CaTi03-Ca(Al 1 / 2 Nbi/2)03 where it was found that the presense of Al decreases 
the Nb dynamical charge [17]. 

The effects of cation ordering on Z* can be well explained via a model in which a 
layered superstructure is treated as a stacking of slabs with different local electronic 
dielectric permittivities Eoo. We fit calculated values to the expressions [2]: 

ZL (Nj) = -^-4,, (3) 

where gij does not depend on the dielectric permittivity. In order to compute 
the dielectric permittivity in the layered structure considered we used a slightly 
modified expression for the permittivity of two sequential capacitors 

£ NTa + Nm (4) 

Nxa/ ZooTa + Kj^Nb/ ^ooNb 

Here ^oo^fe and e^Ta are the electronic dielectric permittivities of pure KNb0 3 and 
KTa03 respectively, iVj (i = Ta, Nb) is the number of Ta- or Nb-planes. This 



expression at A = 1 is explicit for the simple cubic structure. We fixed SooTa at 
the experimental value 5.15 and, correspondingly, \r a ,Ta = 1- Fitting the data in 
Table 1 yields £ooM> = 6.4, gT a ,Ta = 9Ta,Nb = 3.7, gNb,Nb = 9Nb,Ta = 3.5, \m,Nb = 1, 
and \Nb,Ta = ^Ta,Nb = 0.97. Deviation of the nondiagonal element of A from 1, and 
the large values of g, are connected with the difference between the Lorentz field 
in the simple cubic lattice and that in the perovskite structure. The value of the 
electronic dielectric permittivity in KNbC>3 (6.4) is consistent with that obtained 
from first principles computation for a cubic structure [18] (6.6). 

PMN RESULTS 

To further investigate the effects of cation ordering on dynamical charges in 
complex perovskites, Berry phase analyses were performed for various supercells in 
the perovskite based system PbMgi/ 3 Nb 2 / 3 03 (PMN); Results are listed in Table 
4. 



TABLE 4. Dynamical charges Z* ti in PMN: n.d. arc nondiagonal elements where n.d.l cor- 
responds to nondiagonal elements between [0-11] and [110] while n.d. 2 is between [0-11] and 
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As in KTN, the dynamical charges of ions in PMN depend on nn coordination. 
In particular, Z* Ph is sensitive to chemical ordering. These results are consistent 
with Zp fe -computations for some other lead relaxors [20]. It is interesting to notice 
that the largest dynamical charge components for Pb occur in the [111] ordered 
structure. The nondiagonal elements are generally small compared with the differ- 
ences between the diagonal elements, which shows that the basis vectors are close 
to the eigenvectors of the dynamical charge tensor. 

CONCLUSIONS 

We have computed the minimum energy [110] barrier between adjoining [0,0, d] 
Li wells in fully relaxed KLT supercells. Our results are in good agreement with 
values obtained by fitting a mean field expansion to experimental data. 

We find that all ions in the supercell are displaced from their ideal positions 
but among the nn ions, the oxygens closest to Li exhibit the largest sympathetic 



displacements; they are displaced in the opposite direction as Li, and move closer 
to one another. These results support earlier work [11] and computations [12,21]. 

First-principles computation of Li-Li pair excitations in a [K 14 / 16 Li 2 / 16 ]Ta03 su- 
percell strongly supports a previous assumption [15] about the nature of n- 
relaxation in KLT. Our low-energy values for excitations of this point defect are in 
excellent agreement with fits to experimental data. 

Computations of the Born charges in KTN and PMN indicate that dynamical 
charges are sensitive functions of nearest neighbor coordination, which can be ex- 
plained by heterogeneity in the local electronic dielectric permittivity. 
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